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A key open question in the study of life is the origin of biomolecular homochirality: almost every 
life- form on Earth has exclusively levorotary amino acids and dextrorotary sugars. Will the same 
handedness be preferred if life is found elsewhere? We review some of the pertinent literature and 
discuss recent results suggesting that life's homochirality resulted from sequential chiral symme- 
try breaking triggered by environmental events. In one scenario, autocatalytic prebiotic reactions 
undergo stochastic fluctuations due to environmental disturbances. In another, chiral-selective poly- 
merization reaction rates influenced by environmental effects lead to substantial chiral excess even in 
the absence of autocatalysis. Applying these arguments to other potentially life-bearing platforms 
has implications to the search for extraterrestrial life: we predict that a statistically representative 
sampling of extraterrestrial stereochemistry will be racemic (chirally neutral) on average. 



I. INTRODUCTION 

During the past few decades the search for life elsewhere in the universe has risen to the forefront of scientific 
research pQ. This is motivated by recent discoveries of exoplanets, including the discoveries of super-Earths [2] and 
Earth-like planets [3 , opening up the possibility of a potentially large number of habitable planetary platforms beyond 
Earth. In addition, carbon isotopic evidence indicating that life existed on Earth at least as early as 3.5 billion years 
ago (Bya) [H [5] , and the discoveries of extremophilic life forms on Earth [6 , suggest that life can survive and even 
thrive under harsher conditions than previously imagined. In light of such evidence, it is reasonable to conjecture that 
at least simple, unicellular life may be more common in the universe than anticipated. When attempting to answer 
the question of how widespread life is, it is pertinent to examine the only example of abiogenesis known to date: the 
origin of life on Earth. 

One of the most distinctive features of life-the existence of a specific and seemingly universal chiral signature-also 
presents one of the longest standing mysteries in studies of abiogenesis. It is well-known that chiral selectivity plays 
a key role in the biochemistry of living systems: nearly all life on Earth contains exclusively dextrorotary sugars and 
levorotary amino acids. Quite possibly, the development of homochirality was a critical step in the emergence of life. 
Although there are numerous models for the onset of homochirality presented in the literature, none is conclusive: 
the details of chirobiogenesis remain unknown. 

We propose that the environment of early Earth played a crucial role in chirobiogenesis. Environmental effects, 
when strong enough, can destroy any memory of a prior chiral bias, whatever its origin. They can also drive a 
specific chiral choice. Life's chirality is interwoven with early-Earth's environmental history; specifically, with how 
the environment influenced the prebiotic soup that led to first life. The same should be true of life anywhere in the 
cosmos. As remarked in Ref. pQ, the history of life in a planet mirrors the planet's life history. Here, it will be argued 
that the same is true for life's chirality. 

II. LIFE IN THE UNIVERSE 

The rapidity with which life first appeared on Earth is often cited as evidence that life may be common in the 
universe. Observational constraints on the timescale for the origin of life on Earth suggest that Earth-like planets 
older than about 1 Gyr support their own abiogeneses with a probability > 13% at the 95% confidence level [7]. 
Combining this estimate with results demonstrating that Earth-like planets in habitable zones may be a common 
byproduct of star formation [8j[9], the search for extraterrestrial life may yield promising results in the near future. 



* Plenary talk delivered at the Sao Paulo Advanced School of Astrobiology, Sao Paulo, December 2011. 
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However, until extraterrestrial life is discovered, the existence of life on Earth is our only insight into abiogenesis, 
presenting important clues on the likelihood of life elsewhere in the universe. 

Given early Earth's tumultuous environment, it is remarkable that paleontological evidence suggests life may have 
been thriving as early as 3.5 Bya [4j [10]. Perhaps even more surprising are results suggesting that any early life 
may have been killed off as late as 3.8 Bya [TTJ[T2]. These findings indicate that the origin and early diversification 
of life occurred within a window as short as 300 million years. Additional constraints stemming from the half-life 
of prebiotic compounds and the recycling of oceans through hydrothermal vents set an even shorter timescale for 
abiogenesis, estimated to be as brief as 5 million years [T3] . 

The short timescales constraining the origin of life on Earth shed little light on the environmental conditions 
required for abiogenesis. The question of where life began remains open, with few details known about where the first 
prebiotic ingredients may have been synthesized. Potential sites for the origin of life include submarine vents [14], 
primordial beaches [15], and shallow pools and lagoons [16 . An even more debated aspect of the puzzle is whether 
the ingredients for life originated on Earth or were delivered to Earth from outerspace [17] . Evidence supporting the 
latter hypothesis includes the ubiquity of organic chemical ingredients in the interstellar medium [18] and in the solar 
system in carbonaceous chondrites [19 . However, organic compounds have also been shown to be readily synthesized 
under conditions likely to be prevalant on the prebiotic Earth [20H22] . 

Of equal ambiguity is the mechanism of abiogenesis. Although the uniformity of life on Earth suggests that all 
extant organisms descended from a last universal common ancestor (LUC A), we know almost nothing about the abiotic 
ingredients and prebiotic chemistries present on the primitive Earth from which the LUC A evolved [23]. Potential 
mechanisms range from "metabolism-first" models, such as the iron-sulfide world hypothesis of Wachtershauser [24] 
and "membrane-first" lipid- world scenarios as investigated by Deamer and coworkers [25j [26] , to the "peptide-first" 
models proposed by Fox [27, 28 and others [29, 30;, and the popular "genetics-first" hypotheses such as the RNA [3T] 
and pre- RNA [32] world scenarios. In considering any of these models, one must investigate how the characteristic 
properties of life might have arisen-including the emergence of homochirality. 

With such a short window for the origin of life, it is likely that the primordial Earth experienced multiple abiogeneses. 
As has been pointed out by Davies and Lineweaver [33] , if large impacts had frustrated abiogenesis, then as the 
frequency of impacts abated at the end of heavy bombardment, there would have been brief quiescent periods when 
life may have emerged only to be annihilated by the next large impact. Extending this to studies of potential abiogenic 
mechanisms, we must therefore be mindful that in prebiotic Earth reactor pools were submitted to environmental 
disturbances ranging from mild (e.g. tides, evaporating lagoons) to severe (e.g. volcanic eruptions, meteoritic impacts). 
Both kinds of disturbances must have affected the evolution of chirality in early Earth [34-36 . For the remainder of 
this work, we will discuss how such a scenario might have played out. 

III. DECIPHERING THE ORIGIN OF LIFE'S CHIRALITY 

A much-debated question is whether the observed homochirality of biomolecules is a prerequisite for life's emergence 
or if it developed as its consequence [37 3 38 . Adding to the mystery, prebiotically relevant laboratory syntheses yield 
racemic mixtures [39]. This is especially surprising given that statistical fluctuations of reactants will invariably bias 
one enantiomer over the other [40 : every synthesis is ab initio asymmetric. It is thus clear that this asymmetry is 
erased as the reactions unfold. Therefore, a chiral selection process must have occurred at some stage in the origin or 
early evolution of life [41] . 

A common viewpoint is that chiral selection occurred at the molecular level [42J , and that the resultant complexity 
of molecular species led to the eventual emergence of life. This view is based on the argument that life could neither 
exist nor originate without biomolecular asymmetry [38] , and is supported by experiments demonstrating that specific 
conformations of structural entities such as a-helices and /3-sheets can only form from enantiomerically-pure building 
blocks [43j|44] (for an alternative viewpoint see Ref. [35]). Taking this bottom- up approach, we therefore assume that 
the prebiotic conditions necessary for the subsequent development of complex biomolecules had to be chiral. 

A. Modeling Prebiotic Homochirality 

Although it was Louis Pasteur [46], inspired by previous work by Biot, who in the late 1840s was the first to 
recognize that many biomolecules display mirror asymmetry, it was not until the pioneering work of Frank [47 , over 
one-hundred years later, that the first breakthrough in understanding the origin of this asymmetry was presented. 
In this influential work, Frank identified autocatalysis and some form of mutual antagonism as necessary ingredients 
for obtaining biomolecular homochirality from prebiotic precursors. In the ensuing decades, many models exhibiting 
such features have been proposed, each providing its own description of chiral symmetry breaking. 
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The various models presented in the literature range from investigations of simple modifications to Frank's original 
model [48] l5T] , to more recent studies describing the onset of homochirality in crystallization [52] |53], and chiral 
selection during polymerization [54] [55] (see Ref. [56] for a detailed discussion). Among these "Frank" models, one of 
the better known is that of Sandars [54 , which provides a basis for understanding chiral symmetry breaking in a RNA 
world [57] . This model succeeds because it includes both of the necessary features of antagonism and autocatalysis as 
originally proposed by Frank, where the mutual antagonism is provided by enantiomeric cross-inhibition as is observed 
in template-directed polycondensation of polynucleotides [58]. (These terms will be clarified below.) 

As various authors have pointed out [30] [56] [59], when addressing the validity of the Sandars model one must 
consider that the autocatalysis necessary for chiral symmetry breaking in such systems is presently only observed for 
a few non-biological molecules [40j [60] , and would be trying to achieve with even very simple organic molecules [61] . 
Despite this shortcoming, the RNA world hypothesis is still deemed viable by some authors [62 . In addition, the 
Sandars model provides an elegant, and relatively simple, model of chiral symmetry breaking while sharing general 
features in common with other (usually more complicated) models, and as such it has been extensively studied in 
the literature [35] [66] . As we know little about the compositions of the primitive atmospheres and seas [13] and 
even less about prebiotic chemistry [67], it is pertinent to study general features as opposed to details of specific 
models. We thus begin by investigating the Sandars model as the basis of our study with the reasonable expectation 
that the results should be qualitatively similar for other models. We will then briefly describe a recently developed 
polymerization model that obtains a substantial amount of chiral bias without any specific external source: only the 
reaction rates display chiral dependence which may be induced environmentally [63 . 



B. Prebiotic Homochirality as a Critical Phenomenon 

While the details of models describing the onset of prebiotic homochirality mentioned above differ, the qualitative 
features are the same; chiral symmetry breaking occurs due to the introduction of instabilities to the symmetric 
(racemic) state that lead to spontaneous symmetry breaking in physical systems [56] [66]. In other words, the spa- 
tiotemporal dynamics of the reaction network is equivalent to a two-phase system undergoing a symmetry-breaking 
phase transition, where the order parameter is the net chiral asymmetry, A. If we define C and V as the sums of all 
left and right-handed chiral subunits, respectively, then the net chirality may be defined as 

C — V 

Note that the net chirality is symmetric Ao = in the racemic state, and asymmetric A+ y - ^ in the non- 
racemic states. The reaction network is a nonlinear dynamical system with behavior controlled by model-dependent 
parameters, including fidelity of enzymatic reactions [35] [54] El] , ratios of reaction rates [65] , stereoselectivity [68] and 
total mass [30] . 

Just as in other areas of physics, environmental interactions can restore the system to the symmetric (racemic) state, 
even in cases where model parameters are set such that the asymmetric state is stable. In such cases, it is important 
to consider how temperature, or other environmental effects, might work to restore the stability of the racemic state. 
Thinking in this direction, it was Salam [69 who first suggested that there should be a critical temperature, T c , above 
which any net chirality is destroyed. One can think in analogy with a ferromagnet: if heated through the Curie point 
any net magnetization is erased and the system is restored to a symmetric configuration. Here, the net chirality 
plays the role of the net magnetization. While Salam conceded that calculating T c would be challenging using the 
electroweak theory of particle physics (assuming the weak force biases chiral selection [70] [71]), a different route was 
recently taken by Gleiser and Thorarinson [72 . Coupling the reaction network to an external environment modeled 
by a stochastic force, they were able to determine the critical point for homochirality in two and three dimensions. 
We move now to a discussion of their work. 



1. Modeling Spatiotemporal Polymerization 



Although the work of Gleiser and Thorarinson was based on the Sandars model, from the above discussion and the 
work of Gleiser and Walker [65 we expect the results to be quite general. The reaction-network proposed by Sandars 



4 



includes the following polymerization reactions: 

L n + L\ — ? L n+ i, 

Li + L n Di -4 L n +iDi, 

D x +L n D x %D x L n D x , (2) 

supplemented by reactions for D-polymers by interchanging L ^ D (consistent with biochemical usage, we denote 
D-compounds by the letter 11 D" as opposed to the notation set in Sandars' work where such molecules were denoted 
as "R"). A left-handed polymer L n , made of n left-handed monomers, Li, may grow by adding another left-handed 
monomer with a rate fc s , or be inhibited by adding a right-handed monomer Di with a rate fcj. The latter process is 
referred to as enantiomeric cross-inhibition: attachment of a monomer with opposite chirality to one end of a growing 
chain terminates growth on that end of the chain [58] and inhibits any further enzymatic activity. This process is the 
driving force that causes a net asymmetry to develop in this model [65] . 

In addition, the reaction network includes a substrate, S, from which monomers of both chiralities are generated: 
S kc ^ pC ±^ qCD ^ S kc ( pC J^. qCL ^ w h ere p — 4(1 + /) and q = |(1 — /), with / a measure of the enzymatic 
fidelity. Cl(£>) determine the enzymatic enhancement of chirally-pure L(D)-handed monomers, and are assumed to 
depend on the length of the largest polymer in the reactor pool, A/", such that Cx(£>) = Ln(Dn) [54]. Other choices 
are possible, but lead to similar qualitative results [65] . 

Given that the Soai reaction [60] - the most well-known illustration of an autocatalytic network leading to chiral 
purity - features dimers as catalysts [40], we focus on the truncated system for N = 2. We note that it is possible to 
make this truncation while maintaining the essential aspects of the dynamics leading to homochiralization [65 . The 
reaction network is further simplified by assuming that the rate of change of the substrate, [S], and of the dimers, 
[L2] and [-D2]? is much slower than that of the monomers, [L\\ and [Di]. These approximations are known as the 
adiabatic elimination of rapidly adjusting variables [73] and have been shown to produce a reliable approximation to 
the full (n > 2) Sandars model [65] . 

It is convenient to introduce the dimensionless symmetric and asymmetric variables, S = X + Y and A = X — Y , 
where X = [Li](2fc 5 /Qs) 1/2 and Y = [D^fiks /Qs) 1/2 \ respectively [35]. For ki/ks = 1, after a little algebra, the 
reaction network simplifies to 

lo^ = 1"5 2 , (3) 

x dA _ 2JSA 
dt ~ S* + A 2 ' 



where lo = (2k,sQ) 1 ^ 2 has the dimensions of inverse time. S = 1 is a fixed point: the system tends quickly toward 
this value at time-scales of order Iq . 

Substituting S = 1 into eqs. |3| we obtain an effective potential for V(A), 

V(A) = ^--f\n[A 2 + l], (4) 

with fixed points ^4 = 0, ±y/2f — 1. Note that for /< 1/2 an enantiomeric excess is impossible and the only steady 
state is the (stable) symmetric state. In the case / = 1, the potential takes the form of a symmetric double-well, 
where the two fixed asymmetric steady states are homochiral (A = ±1) and represent the global minima. In this 
case, the symmetric state is unstable. 

The form of this potential introduces the possibility of describing chiral symmetry breaking as a phase transition. 
This, in fact, suggests that a proper treatment of the problem should include spatial dependence. To introduce spatial 
dependence to the reaction network, the usual procedure in the phenomenological treatment of phase transitions is 
implemented with the substitution d/dt — ^ d/dt — &V 2 , where k is the diffusion constant [35 . In this coarse-grained 
approach, the number of molecules per unit volume is large enough so that the concentrations vary smoothly in space 
and time. Dimensionless time and space variables are then defined as to = lot, and xo = x^o/k) 1 ^ 2 , respectively. 
For diffusion in water, k — 10 _9 m 2 s _1 , and nominal values ks — 10 _25 cm 3 s _1 and Q = 10 15 cm _3 s _1 , we obtain 
lo = y/2 x 10 _5 s _1 corresponding to t ~ (7 x 10 4 s)to and x ~ (lcm)xo- 

Considering the case where / = 1, for near-racemic initial conditions (\A(0,x,y, z)\ < 10 -4 ), the spatiotemporal 
evolution leads to the formation of left and right-handed percolating chiral domains separated by domain walls, as 
is well-known from systems in the Ising universality class (see Fig. [T]). Surface tension drives the walls until their 
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FIG. 1: Evolution of 2d chiral domains. Red (+1 on the color bar) corresponds to the L-phase and blue (-1 on the color bar) 
corresponds to the D-phase. Time runs from left to right and top to bottom. Top left, the near-racemic initial conditions. Top 
mid and top right, evolution of the two percolating chiral domains separated by a thin domain wall. Bottom left, environmental 
effects break the stability of the domain wall network. Bottom right, subsequent surface-tension driven evolution leads to a 
enantiomerically-pure world. 
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FIG. 2: Time evolution of the spatially- averaged net chirality corresponding to the snapshots shown in Figure^ Stars denote 
times for snapshots and vertical lines mark the beginning and end of stochastic environmental influence. 



average curvature matches approximately the linear dimension of their confining volume. At this point, wall motion 
becomes quite slow, d(A(t))/dt — » 0, where (A(t)) is the spatially-averaged value of the net chiral asymmetry, and 
the domains coexist in near dynamical equilibrium in that the net stresses add to zero (see Fig. [I] top right). The 
time evolution of A(t) is shown in Fig |2l For such model systems, and considering only diffusive processes, it has 
been shown that the presence of a bias from parity- violating weak neutral currents (PV) or most circularly-polarized 
light (CPL) sources [74] (even in the unlikely situation where they could be sustained unperturbed for hundreds of 
millions of years), would not lead to chirally-pure prebiotic conditions [66] . 
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FIG. 3: Average enantiomeric excess versus environmental "temperature" in three dimensions. The error bars denote ensemble 
averages over 20 runs. 



2. Coupling to the Environment: A Critical Point for Homochirality 



Chiral symmetry breaking in the context of this model can be understood in terms of a second-order phase transition, 
where the critical "temperature" is determined by the strength of the coupling between the reaction network and the 
external environment. The external environment is modeled via a generalized spatiotemporal Langevin equation [72] 
by rewriting eqns. [3] as, 

1 -S 2 +w(t,x), 

^(52^42 " l) *), (5) 

where Iq = (2/csQ) 1 / 2 , and w(x.,t) is a dimensionless Gaussian white noise with two-point correlation function 
(w(x! ,t')w(x.,t)) = a 2 5(t' — t)S(x! — x). The parameter a 2 is a measure of the environmental influence. For ex- 
ample, in the mean- field models of phase transitions, it is common to write a 2 = 27/^T, where 7 is the viscosity 
coefficient, ks is Boltzmann's constant, and T is the temperature. Using the dimensionless space and time variables, 
to = lot, and xo = x^o/k) 1 / 2 , introduced above, the noise amplitude scales as a|j — >• Aq 1 (\o/k) d / 2 a 2 , where d is the 
number of spatial dimensions. A crucial point is that even in the case of perfect fidelity, / = 1, where the potential 
supports stable homochiral steady-states, an enantiomeric excess may not develop if a is above a critical value a c . 

As shown in Fig. |3j an Ising phase diagram can be constructed showing that (A) —> for a > a c and chiral 
symmetry is restored. The value of a c has been obtained numerically in two (a 2 = 1.15(/c//o)cm 2 s) and three (a 2 = 
0.65(/c 3 //o) 1//2 cm 3 s) dimensions [72]. Above a c the stochastic forcing due to the external environment overwhelms 
any local excess of L over D within a correlation volume V£ ~ where £ is the correlation length: racemization is 
achieved on large scales and chiral symmetry is restored throughout space. 

In light of these results, it has been shown that within the violent environment of prebiotic Earth, effects from 
sources such as weak neutral currents (which introduce a small tilt in the potential), even if cumulative, would be 
negligible: any accumulated excess could be easily wiped out by an external disturbance [34j[72]. The history of life 
on Earth and on any other planetary platform is inextricably enmeshed with its early environmental history. 



/-1 ^-kV 2 S 



1 dt 

j (d_A 

dt 



kV 2 A 



IV. PUNCTUATED CHIRALITY HU 



The results of the previous section indicate that the environment of early Earth, or other potential prebiotic 
extraterrestrial environments, must have played a crucial role in chirobiogenesis. The chirality of the prebiotic soup 
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FIG. 4: Punctuated Chirality. Impact of environmental effects of varying duration and fixed magnitude (a 2 /a 2 . — 0.96) on the 
evolution of prebiotic chirality in 2d. Short events (last from left), which have little to no effect, should be contrasted with 
longer ones, which can drive the chirality towards purity and/or reverse its trend. (See, e.g. the green line.) 

might have been reset multiple times by significant environmental events such as active volcanism and meteoritic 
bombardment. Under this view, the history of prebiotic chirality is interwoven with the Earth's environmental 
history through a mechanism we call punctuated chirality [34]: life's homochirality resulted from sequential chiral 
symmetry breaking triggered by environmental events. 

Punctuated chirality is an extension of the the punctuated equilibrium hypothesis of Eldredge and Gould [75] to 
prebiotic times. The theory of punctuated equilibrium describes evolutionary processes whereby speciation occurs 
through alternating periods of stasis and intense activity prompted by external influences: the punctuation is the 
geological moment when species arise which may be slow by human standards but is certainly abrupt by planetary 
standards as evidenced by the fossil record. If phyletic gradualism (traditional Darwinian evolution) is like pushing a 
ball up an inclined plane, then punctuated equilibrium is like climbing a staircase [76] . 

It is commonly accepted that molecules undergo selective processes that lead to evolutionary adaptions (see for 
example Refs. [77], [78]). It is therefore natural to extend the punctuated equilibrium hypothesis to the prebiotic 
realm. In this context, the concept of punctuated equilibrium is borrowed with some freedom: the network of chemical 
reactions described in prebiological systems is a non-equilibrium open system capable of exchanging energy with the 
environment. The periods of stasis that develop correspond to steady-states in that even though environmental 
influences may be negligible, chemical reactions are always occurring so as to keep the average concentrations of 
reactants at a near-constant value. 

As an example of punctuated chirality in a prebiotic scenario, one can consider how repeated environmental inter- 
actions influence the evolution of chirality in the context of the model presented in the previous section. In Fig. [4j we 
show several 2d runs where the environmental effects vary in duration, while their magnitude was set at a 2 /a 2 c = 0.96, 
so that the magnitude of the disturbance is just below the critical value found by Gleiser and Thorarinson [72] . Each 
colored line represents a prebiotic scenario, with the same environmental disturbances of different duration occurring 
in sequence. 

In order to investigate the impact of environmental effects on chiral selectivity, the scenarios reflect situations where 
there is no chiral selection, that is, where the two phases coexist in dynamical equilibrium (mathematically, when 
d(A(t))/dt — >• for A(t) 7^ ±1; chemically, in a steady state). We observe that long disturbances can drive the net 
chirality towards purity ((A(t)) —> ±1 for large t). Furthermore, note that subsequent events may erase any previous 
chiral bias, favoring the opposite handedness. In other words, environmental effects of sufficient intensity and duration 
can reset the chiral bias. This is true even if the system evolves toward homochirality prior to any environmental 
event. 

Fig. [5] summarizes the results of a detailed statistical analysis of 100 2d runs that led to initial domain coexistence, 
that is, d(A)/dt ~ (see Fig. [4] for t < 600) [34]. The horizontal axis displays the magnitude of the disturbance in 
units of a 2 . The vertical axis gives the fraction of homochiral worlds, that is, those that after the disturbance obtain 
chiral purity. The colors represent the duration of the event. For a 2 > 0.96a;?, that is, near the critical region, all but 
the shortest events (t < 501q 1 ~ 1.5 months, for the nominal value of Iq = y/2 x 10 _5 s _1 mentioned previously) lead 
to statistically significant chiral biasing. Results in 3d are qualitatively very similar [34] . 
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FIG. 5: Fraction of 2d homochiral worlds after a single environmental event. Colors indicate duration of the event, while the 
horizontal axis labels its magnitude in units of c? c , the critical value for chiral symmetry restoration. The statistical sample 
included 100 runs and the time scales of the events are in units t(/^" 1 ). 



V. CHIRALITY FROM CHIRAL-SELECTIVE REACTION RATES [63] 

Could prebiotic homochirality be achieved without autocatalysis exclusively through chiral-selective reaction rate 
parameters without any other explicit mechanism for chiral bias? In a recent work, we investigated this question, 
focusing on the simplest possible chemical network: no autocatalysis and no enantiomeric cross-inhibition, the trade- 
marks of Frank-like chirobiogenesis models. Instead, we investigated how rare a set of chiral-selective reaction rates 
in a polymerization model needs to be in order to generate a reasonable amount of chiral bias. We quantified our 
results adopting a statistical approach: varying both the mean value and the rms dispersion of the relevant reaction 
rates, we showed that moderate to high levels of chiral excess can be achieved. 

Consider a polymerization reaction network where activated levorotatory and dextrorotatory monomers (L\ and 
D*, respectively), can chain up to generate longer homochiral or heterochiral molecules. (Here, "activated" refers 
to these monomers being highly reactive due to energy input into the reactor pool.) The model reactions include 
deactivation of activated monomers, 

D{ ^ D x (6) 
and the polymerization reactions, 

Dl+D^D j+1 (7) 
Dt+Li^Mn (8) 
D{ + Mij a A M ij+1 (9) 

where ^(D) is the deactivation rate for L(Z))-monomers and a^m) is the polymerization rate for adding activated 
L( J D)-monomers to a growing chain. Here Li and Dj denote homochiral polymers of length i and j respectively, and 
denotes polymers of mixed chirality consisting of i L-monomers and j D-monomers. There is no autocatalysis 
(either explicit or through enzymatic activity of homochiral chains) or enantiomeric cross-inhibition present. We also 
add a source S and a disappearance rate d to the reaction equations so as to model an open system. All reaction 
rates are scaled by the disappearance rate d. 

Instead of assigning ad hoc values for the rates of opposite chirality, we conducted a statistical analysis whereby, for 
different numerical experiments-that is, different realizations of prebiotic scenarios-the values for L and D reactions 
rates were picked from a set allowed to randomly fluctuate about a given mean. In other words, each time we solved 



L\%L X 



L\ + Li — % Li+i 
L\ + Dj a A M X j 
L\ + Mij a A M i+lj 
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the coupled nonlinear ordinary differential equations describing the reactions, we picked a set of random, Gaussian- 
distributed values for the four reactions rates h^o) an d o-l(d)- (For details, see Ref. [S3].) 
We define the enantiomeric excess [ee(i)], that is, the net chirality, as 



ee(t) 



[L\{t) + L x {t) + C(t)} - [Dt(t) + Di(t) + V(t)} 
[Lt(t) + L^t) + C(t)] + [£>!(*) + D x {t) + V{t)} ■ 



(10) 



Note that our definition of the enantiomeric excess implicitly assumes that we are only interested in pure homochiral 
polymers. Here, we have taken the limit of large TV in the polymerization equations and introduced 



(ii) 



for levorotatory, dextrorotatory and mixed polymers, respectively. 

In Figure [6] we show the behavior of the time evolution of enantiomeric excess for several illustrative examples with 
S/d = a = 200 and different choices for the reaction rates a L ^ D )/d = cll(d) an d h L ^ D )/d = Pl(d)- On the left, we 
show results with Pl — Pd = 100 fixed, and varying ratios of oll/old- Shown are the results for ^ = 5.0, cud = 10.0 
(oll/old — \) in blue, oll — 5.0, old — 15.0 [a^jotB — |) in black, and = 5.0, old — 20.0 (oll/old = |) in red. The 
system quickly reaches a steady state with net chiral excess oogg c± 0.24, ee ss ~ 0.38, and ee ss ~ 0.48, respectively 
(corresponding to 24%, 38%, and 48% enantiomeric excess). On the right, we show results holding aL = old = 10.0 
fixed, and varying the ratio of Pl/Pd- Shown are the results for Pl = 75, (3d = 25 (Pl/Pd = 3) in blue, Pl = 100, 
Pd = 25 (Pl/Pd — 4) in black, and Pl — 125, Pd = 25 (Pl/Pd — 5) in red. In this case the net chiral excess 
at steady-state is ee ss ~ 0.38, ee ss ~ 0.47, and ee ss ~ 0.58, respectively (corresponding to 38%, 47%, and 58% 
enantiomeric excess). These test runs show that for differing left and right reaction rates a substantial amount of 
chiral excess may be reached at steady- state. Although for these examples the differences were fairly large, we note 
that if the chiral- selective changes in rates appear in exponential factors (as in the example above for temperature 
dependence), small changes in the parameters may generate fairly large changes in the resulting reaction rates. 
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FIG. 6: Time evolution of enantiomeric excess for illustrative examples of reaction rates with a = 200. Left: Pl — Pd — 100, 
with varying ratios of oil I old- Shown are cxl — 5.0, old — 10.0 (cll/old — \) in blue, cxl — 5.0, old — 15.0 [cxlIold — |) in 
black, and oll — 5.0, old — 20.0 (ollIold — \) in red. Right: aL — old — 10.0, with varying ratios of Pd/Pl- Shown are 
{3 L = 75, Pd = 25 (Pl/Pd = 3) in blue, p L = 100, Pd = 25 (Pl/Pd = 4) in black, and L = 125, Pd = 25 (Pl/Pd = 5) in red 

We have performed a detailed statistical study of the reaction equations with fluctuating values for the L- and 
D- reaction rates oll, &d, Pl and Pd- This involves solving the coupled system of ordinary differential equations N 
times, each with values for the four reaction rates given bya^ = ol + 5l, old = ol + 5d, Pl = P + £l, and P = Pd+^d, 
where the bars denote the mean values and ^l(d) an d £l(d) are Gaussian-distributed random numbers, within a fixed 
rms width set by olo and Pq. 

Our results indicate that the net chirality is overall more sensitive to the amplitude of the chiral- selective variations 
about the mean values of the reaction rates than to the values of the reaction rates themselves. In Figure [7| we explore 
the spread in the distribution of the enantiomeric excess in experimental systems for an ensemble with a = 20, P = 120 
and olo = Po = 0.25. Outliers in the distribution have very high enantiomeric excesses up to as much as 80 — 90%, 
although most systems fall within a range of enantiomeric excesses with \ee\ < 0.25, where the mean of the distribution 
lies. 
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FIG. 7: Distribution of ensemble ee values for N = 5000 with a = 20, P = 120, a = f3 = 0.25, and a = 200 . 

VI. CONCLUSIONS 

During the past decades, several biasing mechanisms have been proposed to explain life's remarkable homochirality. 
Parity violation in weak neutral currents [70], if effective, would provide a universal bias: all amino acids found in the 
universe should be levorotatory. (Although one would still need to explain why sugars are dextrorotatory.) In contrast, 
circularly-polarized UV light [74], if produced in active star-forming regions, would act within a stellar system or, 
at most, within neighboring stellar systems without any uniform bias: in different star-forming regions across the 
galaxy, stellar systems should have stereochemistry with uncorrelated chirality. One of us has recently argued that-if 
relegated to diffusive processes-both mechanisms would probably be ineffective within the time-scales relevant for 
life's emergence on Earth [66 . In any case, the point we are making here is stronger: punctuated chirality would 
render any biasing mechanism ineffective: environmental events have the potential to restore chiral symmetry and 
thus to wash out previous values of chirality locally and, for events of great violence, globally. The same conclusion 
holds for the very different mechanism toward homochirality described in the last section, where no autocatalysis 
or enantiomeric excess is needed to amplify an initial bias; instead, chiral-selective reaction rates alone can lead to 
substantial chiral excess. We noted that this chiral selectivity may also be prompted by environmental effects [63] . 
Both parity violation and circularly-polarized UV light predict that within the same stellar system, each planetary 
platform would have the same chiral bias. In contrast, environmentally-induced chiral selectivity is a local effect, that 
depends on the particular history of the planetary platform harboring life. We thus propose that it is possible to 
distinguish between the three mechanisms through future space missions aimed at studying stereochemistry [79]. If 
chiral bias, as life on Earth, goes through periods of stasis (chemical steady state) punctuated by violent upheavals 
and symmetry restoration, we predict that there would be no chiral correlations even within the same stellar system: 
the same amino acid found, say, in Titan would not necessarily display the same chirality if found on Earth or Mars. 
Of course, only a large enough statistical sample would resolve the issue. 

Biological precursors certainly interacted with the primordial environment and may have had their chirality reset 
multiple times before homochiral life first evolved. In this case, separate domains of molecular assemblies with ran- 
domly set chirality may have reacted in different ways to environmental disturbances. A final, Earth-wide homochiral 
prebiotic chemistry would have been the result of multiple interactions between neighboring chiral domains [35 , 65] [66] 
in an abiotic process that mimics natural selection. 
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